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We generalize a method of control of chaos which uses 
delayed feedback at the period of an unstable orbit to sta- 
bilize that orbit. The generalization consists of substituting 
some portion of the nonlinear dynamical system with a de- 
layed dynamics, rather than using a linear delay function for 
control. A further generalization, in which the control func- 
tion retains memory of all previous periods, allows the region 
of the parameter space over which control can be achieved to 
be extended, but at the price of losing the ability to achieve 
superstability. Nonliner feedback results in a larger basin of 
attraction to the stabilized orbit than linear feedback. For 
a simple test mapping studied (the logistic map) the dimen- 
sion of the system increases from one to two by introducing 
control. We show in the case involving memory, for a particu- 
lar choice of the relationship between the control parameters, 
that the superstable orbit can be recovered without reducing 
the parameter space that can be controlled. This particular 
solution, in addition to having the largest basin of attraction 
of the methods considered, retains the dimension of the un- 
controlled system. 



I. INTRODUCTION 

A number of methods have been proposed for feed- 
back control of chaos Two methods of control 
that stabilize an otherwise unstable periodic orbit have 
received considerable attention recently: 

(a) Ott, Grebogi and Yorke (OGY) ^ introduced a 
method which stabilizes unstable periodic orbits (UPO's) 
found in the chaotic regime via small feedback perturba- 
tions to an accessible parameter. The control pertur- 
bation is given when the orbit crosses a given Poincare 
section, such that the trajectory will be close to the sta- 
ble manifold of the desired UPO. In this method, in the 
limit of zero noise, the orbit of the controlled system 
is identical to the UPO of the uncontrolled system and 
the feedback perturbation vanishes. A drawback for the 
OGY method is that it becomes difficult to apply for 
very fast systems, since it requires computer analysis of 
the system at each crossing of the Poincare section. Also, 
noise can result in occasional bursts where the trajectory 
wanders far from the controlled periodic orbit. 

(b) An alternative method of feedback stabilization of 
UPO's, introduced by Pyragas ||], consists of a continu- 
ous linear feedback applied at each computational time 
step. As in the OGY case, in this method the controlled 
orbit coincides with the UPO of the uncontrolled system 



and the feedback vanishes, for zero noise, when control is 
achieved. The feedback procedure can be applied with- 
out knowing a priori the location of the periodic orbit, for 
a version in which the feedback term contains a delayed 
variable, in which the delay corresponds to the period of 
the UPO. Moreover, it is expected that it can be used for 
fast systems, since no parameters are changed on a fast 
time-scale, and the method does not require a computer 
analysis of the system. For some systems, the method is 
robust even in the presence of considerable noise Q|. A 
disadvantage of Pyragas' method is that it achieves con- 
trol only over a limited range of the parameter space, i.e., 
a given orbit will become eventually unstable in the con- 
trolled system as the parameters are varied more deeply 
into the chaotic regime. The use of delayed feedback 
also increases the dimensionality of the system. Socolar 
et. al. extended the Pyragas method to include, in 
the control term, memory of all the previous states of 
the system, and were thereby able to increase the region 
of the parameter space where control can be achieved. 
The dimensionality of the system also increased in this 
method, to the same extent as in the method of Pyragas. 

The desirable properties of a control system depend on 
the application. Here we consider a system with a fixed 
period UPO which we are attempting to control. The pa- 
rameters that control the coordinates of the UPO may be 
slowly varying compared to the UPO period. The system 
can be considered subject to noise, which may take the 
system coordinates away from the periodic orbit. In this 
general situation we may consider the following proper- 
ties of the control system as desirable: (1) In the neigh- 
borhood of the periodic orbit (assumed stably controlled) 
the actual orbit returns optimally fast to the periodic or- 
bit when perturbed away from it (in the limiting case the 
orbit is superstable). (2) A slow drift in parameters can 
be tracked by the control over the largest possible pa- 
rameter space. (3) In the larger space we (nonlinearly) 
want the stability to be maintained over the widest range 
of uncertainty in either the system parameters or the dy- 
namical variables (basin of attraction) caused by noise. 
(4) In some average sense, we wish to minimize the time 
required to return to the desired solution from the basin 
of attraction (in the nonlinear regime). 

In the next section we examine these criteria for control 
of a fixed point of a mapping corresponding to a periodic 
solution of a continuous system, with a known period. 
We use the well studied logistic map as the test bed for 
the study. Here we introduce modifications to the con- 
trol methods of Pyragas and Socolar, using a nonlinear 
function in the feedback term. As in the OGY, Pyragas 
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and Socolar methods, the stabihzed orbit is identical to 
the UPO of the uncontrolled system, and when control is 
achieved the magnitude of the feedback term vanishes in 
the absence of noise. First we compare the linear control 
of Pyragas to an analogous nonlinear control. Next, we 
compare the case in which memory is introduced into the 
mapping parameters with case in which memory is intro- 
duced into the nonlinear control. In Section III we in- 
troduce a special case of nonlinear control with memory, 
that reduces to a remarkable simple form, whose proper- 
ties are generally better than the other cases examined. 
We show that this latter control can be used efficiently 
to control other periodic points and higher dimensional 
mappings. 



II. COMPARISON OF CONTINUOUS CONTROL 
METHODS 

We start by describing Pyragas' method He con- 
sidered a dynamical system that is governed by ordinary 
differential equations, which are in principle unknown. 
However, some scalar variable y can be measured as a 
system output, and the system also has an input avail- 
able for an external force /. These assumptions can be 
met by the following model, 



(2) 



dy 
dt 



Piy,^) + f{t), 



dx 



Q(2/,x) 



(1) 



where x describes the remaining variables of the dynam- 
ical system which are not available for observation or not 
of interest. The forcing term disturbs only the variable y, 
and it is assumed that the system may be in the chaotic 
regime when the forcing term f{t) is zero. 

Pyragas studied two types of forcing. In the first 
method one determines the UPO yt of the chaotic at- 
tractor from y{t), following well known algorithms 
Then one designs an oscillator that has an orbit equal 
to that of Hi. The forcing term is given in this case by 
f{t) = K[yi{t) — y{t)], where K is an empirically ad- 
justable weight of the perturbation. In the other type 
of forcing considered by Pyragas, the forcing term con- 
tains a delayed term of the variable y, namely, f{t) = 
K[y{t — t) — y{t)], where r is the delay time. If the de- 
lay time coincides with period of the i-th UPO then the 
perturbation f{t) vanishes and y(t) will coincide with 
UPO, as in the first case. However, in this last case, one 
does not need to know the UPO, just its period, nor is 
it necessary to design an external oscillator. [Although 
Pyragas described his method for the situation in which 
one knows only a time series, in all the cases he studied 
the equations that described the system were known.] 
Here we are concerned with the second method, i.e., the 
delayed feedback case. 

Applying Eq. (|l|) for the stabilization of a period-one 
orbit in a one-dimensional mapping F{xn), we have 



The controlled system has dimensionality equal to two 
instead of one for the unperturbed system. The eigen- 
values of Eq. (^) are given by expanding it about the 
equilibrium .x„+i = Xn to obtain. 
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F' + K± .y[F' + K]^ -4:K 



(3) 



where F' = F'{xf) is the derivative of F with respect to 
Xn at the fixed point Xf. 

We illustrate this method of control using the logistic 
map. 



c„+i = F{xn) = 4aa;„(l - a;„). 



(4) 



This map presents a sequence of period doubling bifurca- 
tions as a increases and enters into chaos at a « 0.8925. 
The period-one orbit is stable from a = to a = 3/4. 
The fixed point Xf for the period-one orbit is zero for 
< a < 1/4 and Xf = 1 - l/4a for 1/4 < a < 3/4. 
If a < or a > 1, the attractor is unbounded, diverg- 
ing to infinity. The attractor will also diverge if the 
initial condition xq is not in the interval [0, 1]. The 
period-one orbit loses stability when one of the eigen- 
values has modulus larger than one. For the logistic 
map, for K < 1, an eigenvalue crosses -1, causing the 
appearance of a pitchfork bifurcation. When this occurs, 
Eq. (||) gives F'{xf) = —1 — 2K. Since, for the logistic 
map F' {x f) = 2 — 4a, the bifurcation point a* is 



3 + 2/^ 
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A Hopf bifurcation occurs a,t K = 1, where the eigenval- 
ues cross the unit circle with imaginary values. Beyond 
this value of K there is no stable solution, so the max- 
imum value of a where control can be achieved is given 
by a* = 1.25. 

We compare these results with the use of a nonlinear 
rather than a linear function in the feedback term. Thus, 
the forcing term to stabilize a periodic orbit is given by 
is a nonlinear function G{xn,Xn-i)- Obviously, many 
choices can be made for G, with the constraint that G — 
in the desired UPO. Perhaps the simplest construction is 
G ~ —K[F{xn) — F{x(^n-i)]: with K > 0, since with 
this feedback one does not need to know the equation 
F{xn) that governs the system. We show how this control 
could be applied in the block diagram displayed in Fig. 
This feedback also gave a better performance than other 
nonlinear functions, as discussed below. 

For the period-one orbit, our controlled system can be 
written as 

Xn+l = F{Xn) - K[F{Xn) ~ F{Xn-l)]. (6) 

The eigenvalues for this equation are 
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{l~K)F'±^[{l-K)F'Y+AKr 
-^1,2- 2 ■ 

When a pitchfork bifurcation occurs, for K not very 
large, the most negative eigenvalue is equal —1. From 
Eq. (0) we obtain F'(x/) = -l/(l - 2K). For the logis- 
tic map this gives 



3-4/1: 
A{l-2K)' 



(8) 



From Eq. ^ we see that if X < 0, then a* < 3/4. 
For K > 1/3, the period-one loses its stability not via a 
period doubling bifurcation, but via a Hopf bifurcation. 
For this case, we obtain 



1 + 2K 
AK ' 



(9) 



With Eqs. (|) and ^ we find that the maximum value 
for a where control can be achieved with this method is 
also a* = 1.25, which occurs at K = 1/3.. We com- 
pare a* as a function of K for the nonlinear control in 
Fig. H (solid line) to that using the linear control of the 
Pyragas' method (dashed line). In Fig. ^ we compare 
the values of log2 |A|, with A being the least stable eigen- 
value, for nonlinear control (solid line) with linear feed- 
back (dashed line) and with no control (short-dashed 
line). We note that the transient time (in the linear 
regime), which is proportional to 1/| logj |A||, is smaller 
in the nonlinear control than in Pyragas's method. Also, 
superstable orbit, logj |A| = —00, is preserved with the 
nonlinear control. [The reader might well ask why the 
parameters linearized around the fixed point are not the 
same for linear and nonlinear control. As pointed out 
by a referee, the use of a linear (in x) control parameter 
with K KiAa — 2) brings the results into coincidence. 
However, a dependence of the control on a introduces a 
new complexity into the feedback, which we consider be- 
low]. We note in Fig. ^that the largest range of stable a, 
a < a* , occurs at A' = 1 for the linear control. This has 
significant disadvantages when we consider the nonlinear 
phase space, as we now show. 

We numerically determined the basin of attraction by 
constructing a grid of initial conditions in the Xn-i, Xn 
space and determining which are attracted to the fixed 
point. The effect of noise on the stability is qualitatively 
examined by constructing a noise circle around the fixed 
point, which just touches the basin boundary, and finding 
the radius r of the circle. We illustrate these properties 
in Fig. ^a) for the nonlinear control and in Fig. ^(b) for 
the linear control. In Fig. |^ we show how r varies with a 
for the nonlinear control at AT = 1/3 (solid line) and for 
the linear control along the dotted line of Fig. ^ (dashed 
line). We see that the nonliner control is more robust in 
the presence of noise. 

Although the control procedure we have used is 
straightforward to implement, the reader might well ask 



if a different form of nonlinearity might be better. To 
investigate this possibility we studied the quadratic map 



Xn+i = l-ax^_ 



(10) 



We found that the control term that is most robust in the 
presence of noise and has the smallest transient is given 
hy F = Ka{\xn\'' - \Xn-iY) with z = 2. 

Another method of control of chaos that follows Pyra- 
gas' ideas was introduced by Bielawski et. al. 0. In this 
method the forcing is given to an accessible parameter 
of the system, instead of adding a feedback term to the 
equation. The controlled logistic map in this case is given 

by 



Xn+i = 4(a -I- e„)x„(l - x„). 



(11) 



with e„ = ^{xn — Xn^i)- The method has been gener- 
alized by Socolar et. al. [||, with the controlled logistic 
map given by Eq. (0), but with 
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[Xn - Xn-l) + Re.n-1, 



(12) 



where R < 1. The case i? = reduces to the Bielawski 
control. Socolar et. al. [|| have shown that this form of 
the control parameter is equivalent to including memory 
of all the past states of the system. The dimensionality 
of the new map is also two with the variables a;„ and e„. 
The eigenvalues are 
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F' + R + -i± ^[F' + R + -1Y ~A[F'R + -i] 



(13) 



where 7 — K{Aa — l)/(4a)^. A pitchfork bifurcation oc- 
curs at 

/s: = 8a*^(i?+l)(4a*-3)/(4a*-l). (14) 

A Hopf bifurcation occurs at 

K = 16a*2[l -I- 2i?(2a* - l)]/(4a* - 1). (15) 

The stability boundaries a*{K, R) are given by Eqs. ( |T^ ) 
and (|l5|). By varying R one finds that, in the absence of 
noise, control can be achieved for arbitrary large values of 
a. However, the width of the window in a where control 
can be achieved decreases as a increases. 

Memory can also be included in the form of nonlinear 
control by a generalization of Eq. (^ , 

Xn+l = F{xn) + e„ (16) 

e„+i = -K[F{xn+i) - F{xn)] + Ren, (17) 

When R = this control reduces to the case of nonlinear 
control studied above. The eigenvalues are 
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Ai,: 



(1 - K)F' + R± y/[{l - K)F' + i?]2 + A{K - R)F' 



(18) 



and the stability boundaries for the logistic map are ob- 
tained from 



^ 3{1 + R)-4K 
A{l + R)-8K 

for a pitchfork bifurcation, and 
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1 



2 A{K-R) 



(19) 



(20) 



for a Hopf bifurcation. The maximum value of a where 
control can be achieved in this method is a* — ; 

which occurs at K — -^^q^. In Fig. ^ the stability 
boundaries are shown for the mappings given by Eqs. ( [l^ ) 
and (0) (soHd line) and Eqs. (|ll|) and (|l|) (dashed 
line), with R = 0.5. For both mappings, the addition of 
a memory term (finite R) extends the region in parameter 
space that a can be tracked. There is a distinct difference 
in the results of Fig. ^ for the two methods of control. For 
the Socolar method the range of a that can be stabilized 
becomes small, as K tracks a to large values. In contrast, 
the additive nonlinear control picks out a value of K for 
which the mapping can be controlled for all values of a 
up to a maximum, and thus is not sensitive to parameter 
drift or uncertainty. The range of a that can be controlled 
at the fixed point increases without bound, as i? ^ 1. 

There is a price to pay for having a finite R when the 
variable is subject to noise. For example, at i? = 0.5, 
a = \, K = 0.6428 for the additive nonlinear control 
and K = 1.8333 for the parameter control, which cor- 
respond, respectively, to the diamond and cross symbols 
shown in Fig. ^, we numerically find the basin of attrac- 
tion for the two cases. The result for additive nonlin- 
ear control is shown in Fig. 0(a) and for control in the 
parameter in Fig. 0(b). In both cases the stabilization 
region has been decreased from that without memory, 
but much more so with parameter control for which the 
basin appears to be fractal. The calculation of the noise 
radius here is somewhat subtle, since e does not have a 
clear physical meaning. If we add a noise term 5x to the 
right hand side of Eqs. ( ITTl) and ([l^ ) we note that, when 
the system is at the fixed point, the variation in a:„+i 
will still be 5x. However, the variation in Cn+i will be 
(1 + ^)5x. If the same procedure is applied to Eqs. ( [l^ ) 



and ( |17| ) we find that the variation in e„+i will now be 
(1 — KF' [x f))5x. Therefore the noise 5x is amplified in 
the e variable in both cases, with distinct multiplicative 
factors. To compensate for this, we contract the e coordi- 
nate by the respective factor before calculating the noise 
radius in the basin of attraction of the additive nonlinear 
control and the parameter control. Our results are shown 
in Fig. H for R = 0.5, for nonlinear additive control (solid 



line), calculated at if = 0.6428, and for parameter con- 
trol (dashed line), calculated at the median line of the 
stability boundary shown in Fig. |^. We find that the 
nonlinear additive control is more robust to noise than 
the parameter control. However, for the same value of a 
the nonlinear control with i? > is less robust than the 
case with R = 0, which is shown in Fig. ||. We must also 
consider the effect of varying R. We do this only for the 
case of the additive nonlinear control in the next section. 



III. AN OPTIMAL CONTROL FUNCTION 

Although, the nonlinear control, with a memory fac- 
tor R, has a number of desirable properties, there is the 
drawback that for a given R at the value of K for which 
a* obtains its maximum value, there is no superstable 
orbit (this also occurs in the parameter control). Since 
operation at a superstable orbit is very desirable from 
the perspective of return to the fixed point solution in 
the presence of noise, we look for a relation between K 
and R for which a superstable orbit is recovered by set- 
ting A = in Eq. (fsl). We find two solutions. 



i? = 0, F'{xf) = 0, 

that, for the logistic map, corresponds to a 
is the solution without memory, and 



(21) 



0.5, which 



R = K, K = 



F'{xf) 
F'ixf)-r 



(22) 



which, for the logistic map, corresponds to K 



4a-2 
4a-l ■ 

We call the second solution, Eq. (|22|), an optimized con- 
trol function, as it allows operation with superstabil- 
ity, i.e., with maximum control at the fixed point, for 
< a < oo. We find that this solution has other desir- 
able properties, such as a large basin of attraction, and, 
remarkably, reduces the phase space to a single degree of 
freedom. 

Substituting R = K into Eq. (|l^), and eliminat- 
ing e„ in favor of a;„+i by using Eq. (|l6|), we obtain 
e„+i = -K[F{xn+i) - Xn+i]. Dropping the index by 
1, and substituting for e„ in Eq. ([l^), we obtain a re- 
markably simple form for the mapping equation 



^i=F{x,,)-K[F{x,,)-Xn], 



(23) 



which is valid for control of the period-one orbit in one- 
dimensional maps. At the fixed point, the magnitude of 
the feedback term vanishes, as in the other methods stud- 
ied here. A block diagram of the optimal control scheme 
is shown in Fig. |9[ One expects this method to be easy 
to implement in experiments, since the control term con- 
tains only amplified versions of the input and output of 
the dynamical system and one does not need to know 
F to apply the control. We note that for the particular 
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case of a mapping, the period-one orbit is a fixed point. 
Thus the variable itself can be thought of as a delayed 
signal at the fundamental period of the updated variable. 
This property allows us to use a feedback signal with the 
same index as the mapping function itself. The delayed 
feedback is seen explicitly for control of differential equa- 
tions, as discussed below. The eigenvalue for Eq. ( p^ ) is 
given by 



X = {l^K)F'{xf)+K, 



(24) 



This map loses stability via a pitchfork bifurcation, where 
A = — 1. Consequently, the bifurcation point for the con- 
trolled logistic map is at 



4(1 - K)' 



(25) 



Thus, we see that the parameter region where the period- 
one is stable increases as K increases and tends to infinity 
as K tends to one. 

The superstable orbit, A = 0, is obtained at 



2~K 
4(1 -if)' 



(26) 



where the subscript 's' denotes superstable orbit. Also 
here Us increases with K and goes to infinity as K tends 
to one. In Fig. |l^ we show, as a function of K, the values 
a* where the period-one bifurcates (solid line) and the 
values as of the superstable orbit (dashed line). We plot 
in Fig. |ll| the Liapunov exponent, log2 |A|, as a function 
of a, for K — 0,0.4,0.8. One can see from this figure 
that increasing K increases the range of the parameter a 
around the superstable orbit for which a given transient 
time can be achieved. 

We now calculate the basin of attraction of the con- 
trolled UFO. Since our controlled map is one-dimensional 
this can be found easily. Using Eq. (23) with F given 
by Eq. (Q) the convergence to the UFO will be attained 
when < xo < 1 -|- 4a(^_ji) ■ Substituting for a at the 
superstable orbit from Eq. (|M) we obtain 



0<xo<l 



K 



2- K 



(27) 



The basin of attraction increases with extending from 
to 1 at ii' = to to 2 at if = 1 . Since the fixed point 
isatx/ — 1 — 1/ (4a) , the noise radius around the fixed 
point is 



1 



K 



1 



4a' 4a(l-Jr) 4a 
which, at the superstable orbit, gives 



2 - K' 



(28) 



(29) 



such that Ts varies from 0.5 to 1 as iiT varies from to 
1. Comparing Eq. (|9|) to our previous control parame- 
ters we see that this optimized control maintains stability 
better in the presence of noise. 

Although the superstable orbit is maintained, it is not 
clear what happens to the time constant for return to the 
periodic orbit as, K — > 1, for initial conditions that are 
started far away from the fixed point. To study the effect 
of K on the nonlinear transient we do the following: we 
start the system with 1000 different initial conditions, 
uniformly distributed in the interval [0,1-1- j^f^irj^)- 
Then we verify how many iterations on average are nec- 
essary to bring the orbit within a radius of 10~^ around 
the fixed point. The result of the nonlinear transient as a 
function of K for a = as is shown in Fig. |l^. It increases 
slightly as K increases and goes to infinity at K — 1, 
where Eq. ( p3|) is marginally stable. 

In a more general way, the nonlinear control with mem- 
ory for stabilization of an UFO in a one-dimensional map 
with period m is is given by 



(30) 

F'"(a;„)]+i?e„. (31) 



For the case in which K 
optimized version 



R this control reduces to the 



•^n-^-m (-^n) K[F {Xn) ^n] 



(32) 



Also higher periodic orbits the dimensionality of the con- 
trolled map is still one. The fixed points of the iter- 
ated map are identical to the fixed points of the un- 
controlled equation. We have applied the optimized 
control for a period-two orbit (m = 2) of the logis- 
tic map, in which the fixed points are given hy Xf — 
[4a + 1 ± V(4a-3)(4a+ l)]/8a. The eigenvalue for the 
period-two orbit can be easily calculated and one finds 
that a pitchfork bifurcation from period-two to period- 
four wiU occur when a* = i[l + ^5 + (1 + K)/{1 - K)]. 
The superstable orbit is at as = |;[1 -f \/5 + K/ (1 — A')]. 
The value of a where the bifurcation from period-one to 
period-two occurs is at a = 0.75, which is the same value 
found in the uncontrolled map. Consequently, the region 
of the parameter space where control can be achieved in 
the period- two orbit also grows with increasing K, and 
goes to infinity as K tends to one. We note that the 
period-two orbit is also controllable by the methods con- 
sidered in section ||, but are also "non-optimal" in the 
sense that we have discussed. 

Although the results of our study of controlling a sim- 
ple one-dimensional mapping are suggestive of general 
underlying principles, they are not generic. A generaliza- 
tion of our "optimized" control scheme for a period-one 
orbit can be expressed in the following form 



u„+i = P(u„, v„) f, v„+i = Q(u„, v„). 



(33) 
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where u is a vector of the variables that are available for 
observation and v describes the remaining variables of 
the dynamical system which are not available or not of 
interest. The control term operates only on the u vector 
and is given by f = if[u„ — P(u„,v„)]. We apply this 
more general form to a higher-dimensional mapping. For 
a specific example, we study the Henon map, which is 
given by 

Xn+i = 1 + i/„ - ax^, Vn+1 = bxn. (34) 

In this map, for h = 0.3 (which is the case we con- 
sider here), the period one orbit is stable in the interval 
—0.1225 ^ a < 0.3671. The system enters into chaos 
when a ^ 1.059, and the orbit becomes unbounded for 
a > 1.428. For this map we can use three types of con- 
trol: in both variables, only in the x variable or only in 
the y variable. For the first type of control we have 



Xn+l = 1+ J/n - ax^n - ^[l + IJn " CLX^n ~ ^n], (35) 

y„+i = bxn - K[bxn - Vn]- (36) 
For the second type of control the equations are 

Xn+l = 1 + J/n - axl^ - K[l +yn- CLX^^ - x„], (37) 

y„+i = hxn- (38) 
The third type of control gives 

Xn+l = I + Vn ~ ax'i, (39) 

yn+l ^ bXn - K[bXn ~ Vn]- (40) 



In all these cases the fixed points are the same as in the 
uncontrolled Henon map, and the feedback term vanishes 
when control is achieved. We have found that the con- 
trol does not change the lower boundary of the region 
of stability of the period-one orbit. However, the upper 
boundary increases as K increases for the three types of 
forcing. For example, for K = 0.4 the period-one orbit 
bifurcates at a k, 1.98, 1.74, 0.49 for the the first, second 
and third methods, respectively. Thus, different types 
of control have different regions for which stabilization 
is possible. For the Henon map the largest region of 
control occurs when the x and y variables are controlled 
simultaneously. The largest Liapunov exponent for the 
uncontrolled and for the controlled Henon map is shown 
in Fig. |l^, also for K = 0.4. As we see, no superstable 
orbit exists for any a, for the period-one orbit in the un- 
controlled equation with 6 — 0.3. The feedback terms we 
use to expand the region of stability modifies the location 
of the most stable orbit but do not create a superstable 
orbit when one does not exist for any value of a in the 
uncontrolled equation. 

IV. CONCLUSIONS AND DISCUSSION 

We have generalized a method introduced by Pyragas 
, used to control an otherwise unstable periodic orbit. 



as applied to mappings. The method consists of feeding 
back a delayed signal with the delay equal to the period 
to be controlled, done in such a manner that the position 
of the stabilized orbit in the phase space is not changed. 
The generalization consists of feeding back the nonlinear 
mapping signal, rather than a signal linearized around 
the fixed point. This increases the basin of attraction 
of the controlled signal and thus decreases the sensitiv- 
ity to noise. However, the range of parameters for which 
control can be achieved is limited. An addition to the 
control procedure, introduced by Socolar et. ai Q, is to 
allow memory of all previous periods. This latter proce- 
dure was implemented in the mapping parameter, rather 
than directly into the variables. The method allows an 
arbitrary range of the parameter in the logistic map to 
be tracked, but at the expense of a rapidly decreasing 
the basin of attraction with increasing range of parame- 
ter tracking. A generalization of the nonlinear feedback 
applied to the variables, to include memory, also allows 
arbitrary tracking of the parameter, with a significantly 
improved basin of attraction. All of these above control 
procedures increase the dimensionality of the phase space 
for a one-dimensional map to two. 

For control with memory there are two parameters to 
be chosen, the control parameter K and the memory pa- 
rameter R. For the general case, K , there is no su- 
perstable orbit. However, for nonlinear control with the 
choice R = K the superstable orbit is recovered, and, 
remarkably, the phase space for the controlled logistic 
map is again one-dimensional. Because we recover the 
superstable orbit we call this an "optimized" solution. 

The control methods we have been considering have in 
common with OGY the following properties: (a) the fixed 
points of the controlled map are the same as in the un- 
controlled system; (b) the feedback term vanishes in the 
absence of noise when control is achieved; and (c) one 
does not need to know the mapping equations in order 
to apply the control. Unlike OGY, (d) no computer anal- 
ysis of the system is necessary to apply the control and 
the methods probably can be applied for fast systems; 
and (e) knowledge of the location of the unstable peri- 
odic orbit is not necessary. For the "optimized" control, 
(f) the dimensionality of the controlled equations is the 
same as in the uncontrolled system; (g) the control does 
not destroy the superstable orbit of the uncontrolled sys- 
tem; while simultaneously, (h) control can be achieved in 
a very large region of the parameter space; (i) the basin 
of attraction of the controlled orbit is larger than in the 
other methods; and, consequently (j) the control is more 
robust in the presence of noise. However, (k) to achieve 
control with the parameter values that are deep within 
the unstable region, the nonlinear transient times to re- 
turn to the controlled orbit becomes increasingly long. 

Although we have only considered the application to 
mappings, of the various methods of control, the meth- 
ods are also applicable to continuous systems governed 
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by ordinary differential equations (ODE's). This was 
considered in the original paper by Pyragas who 
applied the linear control to the Roessler, Duffing and 
Lorenz systems. However, unlike a mapping, a sim- 
ple delay makes the dimensionality of the system in- 
finite. We can generalize our nonlinear feedback con- 
trol for the case of ODE's. For the method with- 
out memory the f(t) in Eq. (|l|) is replaced by f(t) = 
-K[P{y{t), x(t)) - P{y(t ~ T),yi{t - r))], where r is the 
period of the UPO. For the nonlinear control with mem- 
ory Eq. (|l|) becomes ^ = P{y{t),x{t)) + e{t), with 
eit) = -if[P(y(t),x(<))-P(2/(t-r),x(t-T))]+i?e(t-T). 
In the case of the "optimized" control we have for the y- 
equation ^ = {1 - K)P{y{t),^{t)) + if^^^. We 
have achieved control of the Roessler system using all of 
these types of feedback. We are currently investigating 
which methods give the best performance with respect 
to the issues that we considered in this paper. 
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FIG. 3. Liapunov exponent, logj |A|, as a function of a 
for the nonlinear control method without memory (solid line) 
with K — 1/3 (the value of K that gives the maximum a*), 
Pyragas' linear control (dashed line) along the dotted line of 
Fig. ^, and for the uncontrolled logistic map (short-dashed 
line) . 



FIG. 4. (a) Basin of attraction and noise circle for (a) the 
nonlinear control without memory for a = 1 and K — 1/3 
and for (b) Pyragas' linear control for a = 1 and K = 0.75. 



FIG. 5. Noise radius r as a function of a for the nonlinear 
control method a,t K — 1/3 (solid line) and for Pyragas' linear 
control along the dotted line of Fig. (dashed line) . 



FIG. 6. Boundary of stability for the period-one orbit 
using our nonlinear control method with memory (solid line) 
and Socolar's control method (dashed line), with R = 0.5 in 
both cases. 
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FIG. 7. Basin of attraction and noise circle for (a) our 
nonlinear control with memory at K = 0.6428, and (b) for So- 
colar's control method a.t K = 1.8333. In both cases, R = 0.5 
and a = 1. These parameters correspond to the diamond 
symbol and cross shown in Fig. ^. 
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FIG. 1. Block diagram for the nonlinear feedback control 
of chaos without memory. 

FIG. 2. Boundary of stability of the period-one orbit for 
the nonlinear control method without memory (solid line) and 
Pyragas' linear control (dashed line). 



FIG. 8. Rescaled noise radius r as a function of a for 
our nonlinear control with memory with K = 0.6428 (that 
is, the value of K that gives the maximum a* in Fig. |^, 
and Socolar's control method along the median line of the 
boundary of stability shown in Fig. ]A In both cases R — 0.5. 



FIG. 9. Block diagram for our optimal control method. 

FIG. 10. a* and as for the logistic map using the optimal 
control method. 

FIG. 11. Liapunov exponent, logj |A|, for K = (solid 
line), K = 0.4 (dashed line), and K = 0.8 (long-dashed line) 
for the logistic map using the optimal control method. 

FIG. 12. Nonlinear transient for the logistic map along 
the line of the superstable orbit, a^, shown in Fig. H. 



FIG. 13. The largest Liapunov exponent, logj |A|, for the 
Henon map without feedback (label 'Henon'), with feedback 
in the x and y variables (label 'xy'), with feedback only in the 
X variable (label 'x'), and with feedback only in the y variable 
(label 'y'). 
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